Vortex images on Bai_a^Ka;Fe2As2 observed directly by the magnetic force microscopy 
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The vortex states on optimally doped Bao.6Ko.4Fe2As2 and underdoped Bao.77Ko.23Fe2As2 single 
crystals are imaged by magnetic force microscopy at various magnetic fields below 100 Oe. Local 
triangular vortex clusters are observed in optimally doped samples. The vortices are more ordered 
than those in Ba(Fei-xCox)2As2, and the calculated pinning force per unit length is about 1 order of 
magnitude weaker than that in optimally Co-doped 122 at the same magnetic field, indicating that 
the Co doping at the Fe sites induces stronger pinning. The proportion of six-neighbored vortices to 
the total amount increases quickly with increasing magnetic field, and the estimated value reaches 
100% at several tesla. Vortex chains are also found in some local regions, which enhance the pinning 
force as well as the critical current density. Lines of vortex chains are observed in underdoped 
samples, and they may have originated from the strong pinning near the twin boundaries arising 
from the structural transition. 

PACS numbers: 74.70.Xa, 74.25.Uv, 74.25.Wx 



INTRODUCTION 



Since the discovery of the iron-based superconductors,^ 
the mechanism of their superconductivity and vortex 
dynamics has attracted considerable interest. Multiple 
electron and hole Fermi pockets, as well as multiple su- 
perconducting gaps, have greatly enriched the physics 
of superconductivity in this new system. Theoreti- 
cally it was suggested that the unique sign-reversal 5- 
wave pairing, namely, s±, could be the main pairing 
symmetry of the iron pnictide superconductors, and the 
nesting between hole and electron pockets is important 
for achieving superconductivity."^ This extended 5- wave 
model results in nodeless superconducting gaps and a 
sign change of the order parameter between the nested 
pockets, which seems to be supported by scanning tun- 
neling microscopy measurement in Fe(Se,Te) samples.^ 
Recent angle-resolved specific heat measurements show 
a fourfold oscillation of the specific heat as a function of 
the in-plane magnetic field direction, which suggests that 
the gap is anisotropic.^ As a type-II superconductor with 
such a multiband and fascinating pairing symmetry, the 
vortex dynamics of pnictides is also attractive. The 122 
family of iron pnictides provides a good opportunity to 
explore the vortex dynamics because of the availability 
of its high-quality single crystals. The parent compound 
BaFe2As2 has both hole and electron pockets with al- 
most balanced charge carriers. Superconductivity can be 
achieved via chemical doping, for example via K substi- 
tution at Ba sites^ and Co substitution at Fe sites. ^ The 
multiband property plays an important role in electric 
transport for both hole- and electron-doped samplesj^iiS 
The magnetization measurements^^^^^ show that the hole 
or electron optimally doped samples both have the second 
magnetization peak effect and a very similar vortex phase 
diagram. For the vortex imaging measurement, almost 



all the measurements were taken on Ba(Fei_ajCoaj)2As2 
or Ba(Fei_ajNia;)2As2 samples. Different detecting meth- 
ods give similar results, i.e., the vortex structure seems 
to be very disordered because of the strong pinning. ^^rJi 
Recently, the scanning tunneling microscopy measure- 
ment on Bao.6Ko.4Fe2As2 shows the ordered vortices as 
well as the Andreev bound states, which seems to be 
very different from Co-doped 122 in which neither the or- 
dered vortex lattice or the in-core Andreev bound states 
were observedi-^^ This is quite natural since the Co doping 
takes place right at the Fe-As planes, while the K dop- 
ing at the Ba sites induces most probably the off-plane 
disorders. Therefore it is very interesting to investigate 
the vortices at low fields on K-doped 122 samples and 
compare them with those in the Co-doped samples. In 
this paper we present the direct imaging of vortices on 
high-quality Bai_a;Ka;Fe2As2 single crystals detected by 
magnetic force microscopy below 100 Oe. The difference 
in vortex structure in both K-doped and Co-doped sam- 
ples is analyzed and discussed in detail. 



II. EXPERIMENTS 

The Bai_a;Ka;Fe2As2 single crystals were grown by 
the self-flux method using FeAs as flux, and the de- 
tailed procedures of synthesizing are similar to previous 



reports. 
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The measurements of x-ray diffraction in- 



dicate a highly c-axis orientation and crystalline quality 
of our samples. The bulk diamagnetic characterizations 
of single crystals were measured by a magnetic property 
measurement system (MPMS, Quantum Design). Mag- 
netic force microscopy (MFM) used in this work is the 
atto-MFM system (attocube) based on the physical prop- 
erties measurement system (PPMS-9, Quantum Design). 
Hard magnetic coating point probes from Nano World 
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were used for all the measurements. The vortex figures 
are made by WSxM software. For every MFM measure- 
ment, the Bai_ajKajFe2As2 single crystal was mounted 
on the sample holder of MFM immediately after it was 
cleaved along the ab plane in air at room temperature. 
The fresh top surface is always fiat and mirror-like for the 
MFM measurements, and usually the measured rough- 
ness on the cleaved surface is less than 1 nm, which ap- 
proaches the measuring precision of the system. Then 
the sample was cooled in a low-pressure helium gas envi- 
ronment. The magnetic property and the MFM measure- 
ments were carried out with the magnetic field perpen- 
dicular to the top surface {ab plane). The magnetic field 
was applied at the temperature above critical tempera- 
ture to obtain the field-cooling process for MFM measure- 
ments. The attocube scanners in the atto-MFM system 
were calibrated by a standard sample to obtain the exact 
scanning parameters at different temperatures. The mag- 
nets of the MPMS and PPMS systems were degaussed 
before the measurements to minimize the residual mag- 
netic field. The first step in the measurement is to find 
a rather fiat place by a tip tapping mode. Then we keep 
a constant distance between the tip and the sample sur- 
face (e.g., 10 nm) and detect the resonance frequency 
change in the presence of the field distribution around 
the vortices. Since the density of vortices changes with 
the magnetic field, we use more scanning pixels to get a 
clearer image at higher fields. 



III. RESULTS 



A. Sample characterization 
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FIG. 1: (Color online) (a) Temperature dependence of volume 
magnetization of the optimally doped Bao.6Ko.4Fe2As2 and 
underdoped Bao.77Ko.23Fe2As2 samples after zero-field cool- 
ing (ZFC) and field cooling (FC) at 20 Oe. The difference 
in the ZFC magnetization of the two samples comes from the 
different demagnetizing effect, (b) Field cooled susceptibili- 
ties versus temperature. Since the susceptibilities are rather 
small, it seems that very few vortices are excluded from the 
sample in the FC process. 



Figure [TJa) shows the temperature dependence of the 
volume magnetization (M) after zero-field-cooling (ZFC) 
and field cooling (FC) processes. Both the optimally 
doped (OP) sample and underdoped (UD) sample used in 
our MFM measurements shows very good superconduct- 
ing transitions. The critical transition temperature of 
the OP sample is 38.8 K (10%Mt=io k) with a transition 
width of 0.6 K, while the value for the UD sample it is 
24.7 K with a transition width of 1.5 K. Bai_a;Ka;Fe2As2 
single crystals are usually very thin, so the demagnetiz- 
ing factor approaches 1.0. That is the reason the ZFC 
magnetization values of these two samples are different. 
The error in the measurement of dimensions, especially 
the thickness, could give error in the calculation of the 
ZFC susceptibility. Figure [Hb) shows the temperature 
dependence of the FC volume susceptibility (xfc), from 
which we can estimate the ratio of vortices excluded from 
the sample. The xfc values of the OP and UD samples 
are only 1% and 6%, which means that a large number 
of vortices are pinned in the samples after field cooling. 
In our MFM measurements we also find that the density 
of vortices is close to that calculated from the magnetic 
field. 



B. Vortex image in OP samples 

In Fig. [2] we show the vortex image on the OP sample 
at different magnetic fields at 2 K by an FC process. The 
maximum field reaches 100 Oe which is almost the limit 
to distinguish the nearest vortices in MFM measurement. 
The calculated number of vortices in this certain range 
was almost the same as the calculated magnetic fiux at 
each field, which is consistent with the FC susceptibility 
mentioned above. The distance between the neighboring 
vortices seems to be very uniform, which is similar to 
the Bitter decoration result on some Ba(Fei_a;Coa^)2As2 
samples, but more ordered than other reports. ^^'^^ Lo- 
cal vortex chains observed in the vortex image are dis- 
cussed in Sec. IIII CI 

To make the figure more clear, we took the coordi- 
nates of all the vortex centers and used the Delaunay 
triangulation to figure out the vortex distribution. The 
result at 100 Oe is shown as an example in Fig. [3l One 
can find that almost half of the vortices are six-nearest- 
neighbored, and more importantly there are distorted 
triangle lattice fragments in some local areas. The self- 
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FIG. 2: (Color online) Vortex image of an OP 
Bao.6Ko.4Fe2As2 single crystal measured by MFM at 2 K and 
different magnetic fields from 10 to 100 Oe. The scanning 
range for each image is 19 /imxl9 /xm. 




FIG. 3: (Color online) Delaunay triangulation of the vortices 
on the OP sample at 100 Oe and 2 K. The red solid circles 
denote 6-neighbored vortices, while the pink empty ones de- 
note the other number of neighbored vortices. To avoid the 
error from the vortices near the scanning edge, only the ones 
with the distance more than 1 /xm from the edges are taken 
into account. The light green blocks show area of distorted 
triangle lattice. 



correlation figures for each vortex image at various fields 
are shown in Figs. Ill^a)-(e). At 10 Oe the irregular loop 
around the center indicates that the nearest distance be- 
tween vortices has a broad distribution, and it means that 
the force among the vortices is very small. The approx- 
imative four fold loop at 20 Oe shows that the vortices 



form a squarelike structure in the local area. With in- 
creasing the magnetic field (or the density of vortices), 
the nearest pattern loop in the self-correlation figure will 
change into a circular shape, which suggests that pairs of 
neighbored vortices have contiguous distances but ran- 
dom orientations. They are not ordered enough to form 
the vortex lattice. Statistics of the distances between two 
nearest vortices by the Delaunay triangulation method 
are shown in Fig. HJ^f). The Gaussian function fits the 
statistic data very well in semilogarithmic scale, and the 
maximum points from the fits are between the values 
calculated from a normal square and a hexagonal vor- 
tex lattice at the same fields. The half- width decreases 
quickly with increasing magnetic field, which means that 
the stronger force at higher fields makes the distance be- 
tween nearest vortices more uniform. It is very difficult to 
get the square pattern in Delaunay triangulation plots, 
so we do the statistic of angles of the Delaunay trian- 
gles as shown in Fig. [5fa). According to the fitting to 
a Gaussian distribution, the square (with characteristic 
angles of 45° and 90°) and hexagonal structure (with 
a characteristic angle of 60°) vortices coexist at 100 Oe. 
Figure [5fb) shows the field dependence of the number ra- 
tio of six neighbored vortices, which behaves linearly in 
a semilog plot. The ratio increases with magnetic field, 
as the vortex system favors a six-neighbored situation at 
high magnetic field. If we do the linear extrapolation on 
the curve to higher fields, the ratio may reach 100 % at a 
field magnitude of several tesla, which means that almost 
all of the vortices have six nearest neighbors at such high 
fields. So it is not strange that the vortices become the 
ordered hexagonal phase in several vortex spacings at a 
magnetic field of 9 T.^^ 

The vortices in the K-doped sample seem to distribute 
more orderly than in the Co-doped sample, which re- 
quires some data to prove. The pinning force per unit 
length and the pinning energy are the parameters char- 
acterizing the vortex pinning. The vortex structure will 
be more disordered if the pinning energy or the pinning 
force is bigger. As a sum of the scalar quantities, the 
pinning energy has a more close relationship with the 
magnitude of the magnetic field than the pinning force, 
so we calculate the pinning force per unit length of the 
vortices to do further analysis. The pinning force per unit 
length on the ith vortex, which is related to the shield- 
ing current at its location from the other vortices, can be 
expressed as 

Here is the magnetic flux quantum, /io is the perme- 
ability of the vacuum, Vij is the space vector between the 
jth and ith vortex, and Ki{rij/X) is the first-order modi- 
fied Bessel function. Figure [6] presents the absolute value 
of pinning force per unit length (/^) distributions for ev- 
ery vortex at 10 and 50 Oe. The statistic result shows 
that most of the values of fi are around 2 x 10~^ N/m at 
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FIG. 4: (Color online) (a)-(e) Self-correlations of vortex im- 
ages at different fields, (f) Statistics of distance between 
the nearest vortices (d) in Delaunay triangulation at various 
fields. The solid line in (f) is the Gaussian fit in semilogarith- 
mic scale. The inset of (f) shows the comparison between the 
peak- value (^center) and the distances expected for a normal 
square (s) and a hexagonal (h) vortex lattices. The errors of 
<icenter are full width at 90% maximum. 



10 Oe, which is 1 order of magnitude weaker than that in 
the Co-doped 122 system at the same magnetic field. 
The penetration depth used here is 0.25 fim at 2 K from 
values of the tunnel diode resonator technology^^ and 
/iSR.^ It should be noted that if we use the larger pen- 
etration depth, i.e., 1.2 /im, as used for the Co-doped 
sample in RefJ^, the calculated value of fi would be- 
come even 1 order of magnitude smaller. So the big dif- 
ference in pinning force in the electron- or hole-doped 122 
system is not from the different chosen values of penetra- 
tion depth. Small-size normal cores are the pinning cen- 
ters in Bao.6Ko.4Fe2As2, which may originate from the 
local doping-induced disorders or some local magnetic 
moments. Although the phase diagrams^^^^^ of super- 
conductivity are similar for both doping sides, the elec- 
tron doping induces impurity by substituting the Fe sites 
with Co, which may be the source of the extra pinning 
centers.^ 

The estimated critical current density jc is propor- 
tional to the pining force per unit length, i.e., jc = //^o- 
Weak pinning in Bao.6Ko.4Fe2As2 means a very small 
critical current density, but it is not consistent with other 
experiment results. One reliable explanation is that the 
K-doped sample has higher Tc, and the pair-breaking 
scattering in the Co-doped sample suppresses the super- 
fluid density. In this way the intrinsic critical current 
densities are very different in K-doped or Co-doped sam- 
ples. It should be noted that the compared pinning force 
mentioned above is the average one. For example, the cal- 
culated average jc = 10^ A/cm^ is very small at 10 Oe. 
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FIG. 5: (Color online) (a) Statistics on the angle values of 
Delaunay triangles. The solid line shows the Gaussian fittings 
at 45, 90, and 60 degrees which are the typical angle of a 
square, or a hexagonal vortex lattice, (b) The proportion of 
the 6-nearest neighbored vortices at different fields in semi-log 
plot. It is obvious that the vortices favor 6- neighbored more 
at a higher fields. The proportion may reach 100% at several 
tesla if we do linear extrapolation in semi- log plot. 



As the penetration depth here (0.25 /im) is much smaller 
than the average distance (1.4 jam) between the neigh- 
boring vortices, the rare vortices have very small interac- 
tions, which may be the reason for the small jc. The peak 
value of jc increases slightly with the magnetic field, as 
shown in Fig.[6fe). At 100 Oe jc at the most frequent po- 
sition is about 1.4 x 10^ A/cm^ and the maximum value 
reaches 1.9 x 10^ A/cm^, which is of the same amplitude 
as the value taken on the magnetization curve. 

The newly cleaved fresh surface is very flat, except 
that some surface steps are formed by the cleaving. In 
Fig. [7] we show the case of two steps with about 10 nm in 
height, and the vortices were pinned by these steps. At 
first glance the pinning of vortices by these steps is sup- 
posed to be induced by the Bean-Livingston pinning. 
When a vortex is close to the parallel mirror surface, an 
attractive interaction is formed between the vortex and 
its image (with opposite sign).^^ In this case the vortices 
should stay at the higher side of the stage. This kind of 
pinning can only happen when the step is high enough 
leading to a large mirror area parallel to the vortices. 
This can readily explain why the vortices along the up- 
per step locate on the high stage and keep some distance 
away from the step, not on top of the clear-cut line. How- 
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FIG. 6: (Color online) Two-dimensional color mapped pin- 
ning force per unit length of every vortex at 10 Oe (a) and 
50 Oe (c), while the distributions of the pinning force per 
unit length are shown in the histograms of (b) and (d), re- 
spectively. The highest frequency value of the pinning force 
per unit length at 10 Oe is about 2xl0~^ N/m, which is 1 or- 
der of magnitude smaller than the value in the Co-doped 122 
system at the same field.— (e) Pinning force per unit length 
and the related critical current density distributions at each 
magnetic field. 




FIG. 7: (Color online) Vortex image (a) and its differential 
image (b) in the presence of two steps on the surface at 10 Oe 
and 4.2 K. The insert in (a) shows the height landscape along 
the white line. It is clear that the heights of the two steps 
are about 10 nm, and the image dimension is 19 /imxl9 /xm. 
The vortices near the steps were pinned along the edge of the 
steps. 
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FIG. 8: (Color online) Vortex image (a) and its self- 
correlation map (b) of an underdoped Bao.77Ko.23Fe2As2 
sample at 10 Oe and 2 K with a scanning range up to 
22 /imx22 fim. The white line shows the main distance of 
the chains, which is about 3.2 /xm in each small division. 



ever, the vortices along the bottom step seem to locate in 
both sides of the line. This may suggest that twin bound- 
aries are induced near the step, which may have weaker 
superconductivity and thus construct a strong pinning 
well. This is similar to the recent report of the vortex 
state near the twin boundaries. 



C. Vortex image in UD sample 

As mentioned in Part IIIIB) there are some vortex 
chains in optimally doped sample as shown in Fig. [2l 
while the pinning force per unit length is also large on 
these vortices, as shown in Figs.[6fa),(c). There are some 
reports on the vortex chain state in the presence of a 
tilted magnetic field from the c-axis in cuprates^^^^^. It 
was shown later that these vortex chains are formed by 
the vortex pancakes which are dragged by the underneath 
Josephson vortex. Therefore to observe the vortex chains 



in Bi-2212, we need a significant misalignment, i.e. more 
than 45°, between the direction of the magnetic field and 
the c-axis of the sample, while in our experiment the mis- 
alignment is smaller than 5°. Additionally, no evidence 
for either pancake vortices or the Josephson vortices was 
observed in the K-doped 122. 

To investigate the vortex chains, we did the further 
experiment on the underdoped samples. Fig. [8] shows 
the vortex image of the underdoped Bao.77Ko.23Fe2As2. 
Clearly, there are vortex chains along the same direc- 
tion, and the distances between the chains are sev- 
eral micrometers. The vortex chains in underdoped 
sample are very similar to the ones found in twinned 
YBasCusOr-^^^'^^ or ErNisBsC^^. The iron pnictides 
have the spin-density-wave and the structural transition 
in underdoped samples^^^^^. This structural transition 
from orthorhombic to tetragonal causes the twin bound- 
aries which are parallel each other and have several mi- 



6 



crometers in distance^. The self-correlation result of 
these vortex images shows that the averaged space be- 
tween the two vortex chains is about 3.2 /im, which is 
consistent with the space between the domain walls. No 
doubt, these twin boundaries enhance the critical cur- 
rent density as counted from the gradient of the vortex 
density. ^^'^^ If we come back to the case in optimally 
doped sample, some strong pinning centers also show up 
near the vortex chains, as shown in Fig. [6l Therefore, 
even in the optimally doped sample, we also have the 
twin boundaries as the strong pinning centers. It is these 
strong and large scale pinning centers that enhance the 
critical current density greatly in the weak field region, 
leading to a sharp magnetization peak near zero field. 
When the magnetic field is increased to a high value, 
more vortices fill into the area between the "network" of 
these twin boundaries. Therefore it is quite natural to ob- 
serve some complex structure of magnetization hysteresis 
loops which exhibit multiple magnetization peaks. Our 
MFM data here give a direct visualization of this kind of 
picture for vortex pinning in the iron pnictide supercon- 
ductors. 



IV. CONCLUSIONS 

We present the vortex images on optimally doped 
Bao.6Ko.4Fe2As2 and underdoped Bao.77Ko.23Fe2As2 sin- 
gle crystals with magnetic fields below 100 Oe. The 
vortices are very diluted and widespread when the mag- 
netic field is several oersteds, but they get crowded with 



a similar distance at higher magnetic field (higher than 
20 Oe). Some vortex chains are observed together with 
a roughly random distribution (with short-range hexag- 
onal order) of vortices between them. The calculated 
pinning force per unit length seems much smaller than 
that in the Co-doped 122 system at the same field, indi- 
cating that the pinning at Fe sites yields stronger pinning 
and vortex disorders. The vortex system becomes more 
ordered and favors a six-neighbored structure at higher 
magnetic field. We find some surface steps as the pinning 
centers but they may not act as the Bean-Livingston mir- 
ror pinning. A vortex chain state is also observed in the 
underdoped sample and is ascribed to the pinning by the 
twin boundaries generated by the structural distortion of 
the orthorhombic state. We observed a cooperative pin- 
ning induced by the large-scale twin boundaries and the 
weak local disorders, which may be a common picture to 
describe the vortex dynamics in iron pnictide supercon- 
ductors. 
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